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Several hitherto unexplained features of thermal conductivity in group IV materials, such as the change in 
the slope as a function of sample size for pure vs. alloyed samples and the fast decay in thermal conductivity 
for low impurity concentration, are described in terms of a transition from a collective to kinetic regime in 
phonon transport. We show that thermal transport in pure bulk silicon samples is mainly collective, and that 
impurity/alloy and boundary scattering are responsible for the destruction of this regime with an associated 
strong reduction in thermal conductivity, leaving kinetic transport as the only one allowed when those resistive 
scattering mechanisms are dominant. 
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I. INTRODUCTION 


It was already recognized by Peierls that normal scatter¬ 
ing is the key point to understanding thermal conductivity [1]. 
However, the preference for a kinetic point of view and the in¬ 
sufficient data on micro and nanoscale samples in those years 
led to the success of Callaway’s proposal [2], and therefore 
the need to treat normal and resistive processes on a different 
footing was never properly acknowledged. It has only been 
in recent years, with the appearance of thermal transport mea¬ 
surements at the nanoscale [3], that the need for a deeper un¬ 
derstanding has become evident, prompting some authors to 
work with modified expressions for the relaxation times, con- 
hnement effects or including new scattering mechanisms in 
the Callaway model [4, 5]. 

Solutions of the Boltzmann equation with ab initio scatter¬ 
ing rates with the proper inclusion of the role of normal (N) 
scattering have achieved in the last years a great level of ac¬ 
curacy, and the properties of large number of materials have 
been predicted without the need of any fitting parameter [6, 7]. 
However, ab initio results are much less amenable to analy¬ 
sis and the computational cost may become unacceptable for 
samples with low dimensionality, showing alloy scattering, 
and/or strong surface effects. To overcome that, analytic ex¬ 
pressions have been proposed from the ab initio calculations 
for the relaxation times in order to be used in simple approxi¬ 
mate models [8]. 

In recent papers a different approach, the Kinetic Collective 
Model (KCM), has been proposed [9, 10]. The model is based 
on a combination of a Kinetic and a Collective contribution, 
weighted by the relative importance of normal vs. resistive 
processes. The main differences with a pure kinetic model are 
the introduction of a hydrodynamic (collective) flow which 
has a unique value for the relaxation time for all the modes 
and uses a form factor to include boundary effects. 

Recent experimental works have demonstrated that pure ki¬ 
netic models are not enough to understand thermal conduc¬ 
tivity at short scales and times [11, 12]. In these works, 
collective flow has been used to explain the origin of the 
non-monotonous dependence of the thermal boundary resis¬ 
tance as a function of the size of the heater arising from ul¬ 


trafast laser heating experiments [11]. Also theoretically, 
collective transport has been successfully used to understand 
first-principles results on graphene thermal transport [13, 14], 
showing that Poiseuille flow can be behind the high thermal 
conductivity of this material and thus converging with the 
KCM. Collective models will be necessary in next years in 
order to analise these experiments and simulations. 

This Kinetic Collective Model has been able to predict ther¬ 
mal conductivity of natural and isotope rich silicon bulk, and 
micro and nanoscaled samples by fitting only the natural bulk 
sample [9]. In a following work [10], the calculations were 
extended to all Group IV materials, signiflcatively extending 
the predictability range of previous analytical models. A no¬ 
ticeable point is that the expressions for the relaxation times 
that seem to work better with the KCM are those obtained 
in the works by Herring and Klemens [15, 16], pointing to 
the possibility that the problems appearing in the last decades 
when predicting thermal conductivity may not be related to 
the relaxation times, but rather to the insuficient attention to 
this collective regime. 

In this Letter we go a step forward in the KCM and predict 
the thermal conductivity in Si/Ge bulk and nanowire alloys 
from the same natural silicon fitting already used in Ref. 10. 
From our predictions we demonstrate that the difficulty of pre¬ 
dicting the thermal conductivity in this kind of samples is re¬ 
lated to the change from a collective to a kinetic dominated 
flow. The fact that such a variety of materials, isotope frac¬ 
tions, alloy compositions, and length scales is correctly de¬ 
scribed by the KCM is a very strong indicator of its validity 
as an analysis, prediction and optimization tool. 

From the KCM, the lattice thermal conductivity 

K = Kk(l - S) -f Kc7*’(Leff)E (1) 


is obtained as the sum of the kinetic Kk and the collective Kc 
contributions weighted by a switching factor 


E = 


1 


(l + (rA)/(rfi)) 


G [0,1] 


( 2 ) 


This factor determines the actual fraction of heat carried 
in each regime depending on the strength of normal (r^v) 
and resistive (Tp) scattering times. Also, the form factor 
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F{Les) G [0,1] modulates the hydrodynamic reduction of 
the flux on the collective term due to boundary effects, like in 
a Poiseuille flow, where Lgs is the effective size of the sam¬ 
ple. Expressions for /tk and Kc, E and F{Leg) are given else¬ 
where [9, 10]. 

In order to compute the properties of Sii-ajGea, alloys, we 
need the phonon dispersion relations and the relaxation times. 
To obtain the dispersion relations for a stochiometry x, we 
perform lattice dynamics calculations within the Virtual Crys¬ 
tal Approximation (VCA) in an ab initio framework [17]. Cal¬ 
culations were done with the Quantum ESPRESSO pack¬ 
age [18], which implements de Density Eunctional Theory 
(DET) [19, 20], under the Local Density Approximation in 
the parametrization of Perdew and Zunger [21]. Core elec¬ 
trons were accounted for with norm-conserving pseudopoten¬ 
tials of the von Barth-Car type [22, 23], and plane waves were 
cut off at an energy of 60 Hartree. For each composition x, the 
lattice parameter was adjusted until the pressure was less than 
0.1 kbar. Solution of the ensuing dynamical matrix provided 
the sought dispersion relations and derived quantities. 

In this kind of alloys four different relaxation times should 
be considered; three of them are significatively contributing to 
the total resistivity, and they should be combined using Math- 
iessen’s rule 



TR Tj Tb tjj 


(3) 


where tb is the boundary term, tjj the umklapp contribution 
and T/ the impurity term. The fourth element is the normal 
scattering term tat. All these terms may depend on x and the 
phonon dispersion. 

For boundary scattering we use 

= Eeff/fx (4) 

where L^g is the characteristic length of the sample and Vx is 
the group velocity for stoichiometry x. 

For Umklapp and Normal scattering we use the same ex¬ 
pressions as in [9]: 


'^N,x 


1 ^ 1 
Bn,xT ^ - expi-ST/Qo)] 


(5) 


and 


^ 1 , _ expjQu/T) _ 

Bu,x 0 o^T[l - exp(-3r/0c)] ' ^ ^ 

where 0£) is the Debye temperature, Qjj the Umklapp freez¬ 
ing temperature , and Bjj^j^ x the composition dependent pa¬ 
rameters, which can be calculated from those of pure silicon 
Busi with [24, 25] 


B{u/N)x — f(U/N)xB[u/N)Si 
B[u/N)x = f[u/N)xB[u/N)Si 

and using only material properties 




(7) 

(8) 

( 9 ) 



FIG. 1: Thermal conductivity in terms of Ge fraction x. Green filled line is the kinetic 
and red line the collective fraction of the thermal conductivity. The top of the 
green line is the total thermal conductivity. Crosses are the experimental 
thermal conductivity obtained from [26, 27] 


Following this procedure, the only fitting values are those 
of pure Silicon: Busi = 2.8 • 10“^® s^K“^ and B^Si = 
3.9 ■ 10-23 s-iR-i, = 7-103 s-iR-i and B'^g^ = 
7.5 ■ 10® s-^K”^. The rest of the values are calculated from 
the dispersion relations. All the independent parameter values 
for the alloy can be computed from linear interpolation, except 
for the phonon velocity in Eqs. (9), which requires Vx = ((1 — 
x)Vs‘^ +XVqI)-^/‘^. 

The alloy relaxation time needs a more detailed discussion. 
In single specie crystals the mass defect term describes the 
variability in isotopic abundance but in alloys it should also 
account for the variability in the force and lattice constants. 
Thus, including these three factors, the total impurity relax¬ 
ation time can be expressed as [16]: 

- = jVS^uj^Dx ( 10 ) 

Tj b 

where as in Ref. 28, instead of the Debye approximation, we 
use the DOS obtained from the full dispersion relations at sto¬ 
ichiometry X, Dx, and the variance term from Ref. 16 


^ + gr„2 + ^—Efl, 


where 


Ta — ^ ^ 

i 



( 11 ) 


( 12 ) 


is the coefficient of variance of (M) mass, (v^) squared ve¬ 
locity or (i?) impurity radius, being ai the value for the iso¬ 
tope/species i and ax the averaged value over all the atoms, 
7 = 1.7 is the experimental Griineisen parameter and Q = 4 
is a factor depending on the geometry of the impurity (substi¬ 
tutional in this case). 

As expected, the second and third terms reduce to zero for 
pure silicon or germanium and the isotopic mass defect corre¬ 
sponding to the first term is the only one remaining. 
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FIG. 2: Theimal conductivity in terms of size Left. Upper plot correspond to pure Si 
and lower plot to Sio. s Geo, 2 . The green filled line is the kinetic fraction and 
the red line the collective fraction of the thermal conductivity. The black line is 
the total thermal conductivity. 


Notice the simplicity in the relaxation times employed (4)- 
(11), where only compositionally averaged values for the pa¬ 
rameters are used. 

Figure 1 shows the good prediction given by the model for 
bulk alloys (7 mm wide rods) at 300 K. The black line over 
the filled area shows the total thermal conductivity, and the 
green and red areas show the kinetic and collective fraction of 
the total transport, respectively. Note that collective transport 
is only important for very pure materials (Si or Ge), with im¬ 
purity fractions of the order of 0.5% already destroying most 
of the collective regime. The prediction given by the current 
proposal in this region is in remarkable agreement with ex¬ 
perimental results from [27], showing that in order to have a 
model able to predict at the same time the thermal conductiv¬ 
ity for pure and alloyed materials, proper description of the 
evolution of the collective contribution is needed, which is 
achieved in the KCM by the inclusion of S. This might be the 
reason why pure kinetic expressions need an extra parameter 
in the alloy term [29]. 

In Fig. 2 we represent the kinetic and collective contribu¬ 
tions for two different situations, for pure silicon at the top 
and Sio. 8 Geo .2 at the bottom. In pure Si, it can be observed 


that in the region where the thermal conductivity is not af¬ 
fected by size effects, the kinetic/collective fraction does not 
change significatively. When boundary effects are noticeable, 
the reduction in the collective flux is more important and the 
ratio kinetic/collective raises. This is because boundary scat¬ 
tering is resistive and its appearance reduces the weight of 
the collective term. We notice that the disappearance of the 
collective contribution is the responsible for a much stronger 
decrease in thermal conductivity in the Lett =1 /rm-10 nm 
range compared to the micron and above scale. [37]. Present 
results agree with recent works pointing on the direction that 
long mean free path phonons are important to predict ther¬ 
mal conductivity in alloy samples [30-32]. For the alloy sam¬ 
ple, in the bottom figure, the collective term has been com¬ 
pletely eliminated due to the impurity scattering and boundary 
scattering does not change the regime. Also, the functional 
dependence with sample size is completely different in both 
samples. The alloy shows the divergence expected for 
kinetic samples dominated by boundary scattering [29] while 
the change in the slope in the pure Si case seems to be behind 
the destruction of the collective term regime, leaving a linear 
dependence with size at low temperatures. High-frequency 
phonons in the collective regime drag low frequency phonons, 
and the increase with sample size is much reduced at high tem¬ 
peratures [33]. 

To show the predictive capability of the KCM in nanome¬ 
ter alloyed samples, in Figure 3 we have plotted the thermal 
conductivity for bulk silicon and germanium, as well as a se¬ 
lection of Si/Ge alloy nanowires. All the curves, showing a 
remarkable agreement with the experimental points, are cal¬ 
culated from the same single fit to natural bulk silicon, in¬ 
dicating that the change in the collective/kinetic ratio seems 
to be the key in predicting this kind of attribute. For pure 
Si nanowires [3], alternative interpretations based on a purely 
diffusive viewpoint have been proposed[34, 35]; nevertheless 
the KCM has been able to adjust the results by Li et al. [10] 
as well as the SiGe alloy nanowires in Fig. 3. While sur¬ 
face roughness in nanowires may have an important impact 
on thermal conductivity [34, 35], we haven’t included it in our 
analysis because the wires in Ref. 36 had sub-nm roughness. 

In conclusion, thermal conductivity values derived from the 
Kinetic Collective Model for pure Si, Ge and alloy Si/Ge sam¬ 
ples ranging from bulk to nanowires for a wide range of tem¬ 
peratures show excellent agreement with experimental data, 
based on a single fit to natural bulk silicon. No further ad¬ 
justable parameters are needed, showing that the difficulty in 
previous attempts to fit wide ranges of temperature, sizes and 
compositions seems to be related with the transition from the 
collective to the kinetic regime experienced by a sample with 
increasing resistive scattering. This, together with the fact that 
the expressions for the relaxations times agree with theoretical 
predictions without any modification, gives plausibility to the 
model. Boundaries and species variation reduce the collec¬ 
tive transport in alloys, with alloy scattering a more effective 
mechanism than boundary scattering. Alloy concentrations as 
little as 0.5% or system sizes of the order of 10 nm will de¬ 
stroy the collective contribution to the thermal conductivity. 
This insight could be very relevant when addressing phenom- 
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FIG. 3: Thermal conductivity in terms of Temperature for samples with Ge content and 
size. Curves are computed according to Eqs. (l)-(l 1). Experimental data are 
obtained from [36] 


ena such as phonon drag or dopant effects in semiconductor 
thermoelectrics. 
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